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Abstract: We study the reliability of the Two-Step moduli stabilization in the type-IIB 
Large Volume Scenarios with matter and gauge interactions. The general analysis is based 
on a family of A/" = 1 Supergravity models with a factorizable Kahler invariant function, 
where the decoupling between two sets of fields without a mass hierarchy is easily un- 
derstood. For the Large Volume Scenario particular analyses are performed for explicit 
models, one of such developed for the first time here, finding that the simplified version, 
where the Dilaton and Complex structure moduli are regarded as frozen by a previous 
stabilization, is a reliable supersymmetric description whenever the neglected fields stand 
at their leading F-flatness conditions and be neutral. The terms missed by the simplified 
approach are either suppressed by powers of the Calabi-Yau volume, or are higher order 
operators in the matter fields, and then irrelevant for the moduli stabilization procedure. 
Although the power of the volume suppressing such corrections depends on the particular 
model, up to the mass level it is independent of the modular weight for the matter fields. 
This at least for the models studied here but we give arguments to expect the same in 
general. These claims are checked through numerical examples. 

We discuss how the factorizable models present a context where despite the lack of a hierar- 
chy with the supersymmetry breaking scale, the effective theory still has a supersymmetric 
description. This can be understood from the fact that it is possible to find vanishing 
solution for the auxiliary components of the fields being integrated out, independently of 
the remaining dynamics. 

Our results settle down the question on the reliability of the way the Dilaton and Complex 
structure are treated in type-IIB compactifications with large compact manifold volumina. 
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Introduction 





The great progress that string theory phenomenology has achieved in recent years is in 
contrast with the several approximations these studies have always implicit, being one of 
the most drastic to neglect a large subset of fields. The first candidates to be neglected are 
the string excitations and Kaluza-Klein (KK) modes which are always very heavy fields, 
with masses dictated by the string and compactification scales. Then, decoupling [1] and 
symmetry arguments can be used to argue for an "effective" theory, although these have 
not been properly integrated out [2]. At this level the truncation is fine since these fields 
and their couplings are easy to be identified. Nevertheless, the Supergravity (SUGRA) 
theory describing the remaining light modes is still very complicated such that an explicit 
and precise study continues to be extremely difficult. A further truncation is then compul- 
sory but now over fields whose masses and couplings are highly sensitive to the point we 
stand in the field space, so that this procedure starts to be less clear. 
This is precisely the philosophy adopted in the stabilization of moduli fields in string 
compactifications. In this context the moduli to be fixed by tree level effects in the su- 
perpotential, fiux induced terms [3] which are quantized and therefore naturally of order 
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one in Planck units, are regarded as completely frozen. Then the rest of the fields are 
stabilized using non-perturbative effects, that are naturally suppressed. As presented this 
Two-Step procedure seems alright due to the hierarchy between the sources of the dynam- 
ics stabilizing each sector; however, the moduli couplings and even their masses are not 
as easy to spot as in the case of the string and KK modes, then a more careful study is 
required. In recent years several studies has tackled this issue from different perspectives 
(see for example [4-16]). Particularly in [10, 12] things were settled down for models where 
the perturbative part of the superpotential takes a tiny value at the vacuum. There it was 
pointed out that this tuning on the superpotential, firstly motivated in the seminal work 
of Kachru-Kallosh-Linde-Trivedi (KKLT) [17] as such to obtain good phenomenological 
Vacuum Expectation Values (VEV) for the moduli, is also a necessary requirement for the 
reliability of the Two-Step procedure, as it ensures a mass hierarchy between the two sets 
of fields. These results apply to the several sequels of KKLT where the VEV of the super- 
potential is quite small, and the fields to be fixed by the fiuxes are completely disregarded 
focusing only on the second step where the remaining fields are stabilized. 
An important point here is the fact that the theory used in the second step is still a A/" = 1 
SUGRA theory, so the first step is supposed to proceed in a supersymmetric (SUSY) way. 
Prom this point of view the tuning on the superpotential can be understood using the fact 
that the SUSY breaking scale in these KKLT scenarios is dictated mainly by this VEV, 
so that being small the heavy fields decouple from SUSY breaking effects and the effective 
theory is approximately SUSY [13]. 

Pollowing this logic, things seems more subtle for models where the neglected fields acquire 
masses that are of the same order of the masses of the fields to be stabilized and the SUSY 
breaking scale. In this class of models fall the Large Volume Scenarios (LVS) [18] where 
the masses and scales are mainly ruled by the size of the compact manifold developing 
exponential sized volume. 

The LVS has become one of the paradigms of string moduli stabilization being widely 
studied (see for instance [19-23] and more recent progress in [24-29]). The main appealing 
feature of this scenario is the fact that it is not necessary to claim for a tuned VEV for 
the superpotential in order to obtain a low SUSY breaking scale. However, as usual in 
any moduli stabilization scenario in type-IIB orientifold compactification, the Dilaton and 
Complex structure moduli are regarded as frozen although their masses lie on the same 
scale off the remaining field masses, as can be understood from the SUGRA contributions 
to the masses, coming like ~ (e^|VF|^), with K the Kahler potential and W the super- 
potential. Still, from the structure of the scalar potential, in the large volume limit, it is 
possible to argue for a decoupling [18] and systematic treatments on this issue have been 
worked out [8, 10]. The idea in this cases is the realization of a weak coupling between the 
two sectors so that the stabilization of one side be almost insensitive to the other one. 
In absence of gauge interactions such a decoupling is achieved if the system with two sec- 
tors, identified by H , the "heavy" fields, and L, the remaining ones in the effective theory, 
is described by a Kahler invariant function, G = K + Mp log (| Vl^|^/Mp) , Mp the reduced 
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Planck mass, with a factorizable form, i.e., 

G{H, H, L, L) = GniH, H) + Gl{L, I) . (1.1) 

This scenario was proposed by Binetruy et al. [30] and extensively studied by Achucarro 
et al. in several works [6, 8, 9, 14, 16]. In [8, 10] it was pointed out that a Two-Step pro- 
cedure in the LVS scenario, for the pure moduli case, can be validated through the same 
generic arguments since all mixing terms in G between the Kahler sector and the Dilaton 
and/or the Complex structure moduli are suppressed either by powers of the volume or by 
non-perturbative effects. Despite few comments in [9] these studies where restricted to the 
case where no gauge interactions nor matter fields were involved. Remains, then, and is 
the aim of this paper, a systematic study for these more realistic scenarios. 
We will follow the same philosophy of [10, 12] by performing a proper integration of the 
H fields and compare the resulting theory with the simple naive theory were the H fields 
are simply frozen. Recently Brizi et al. studied the way fields can be integrated out in 
a SUSY fashion for SUGRA theories [13]. One of their results is that as far as there is 
a hierarchy between the mass of the fields integrated out and the SUSY breaking scale, 
one can neglect these breaking effects and integrate the fields in a SUSY way, so that 
the effective theory, at first approximation in a derivative expansion, is still a SUGRA 
theory. They, moreover, found that the integration of chiral multiplets proceeds, again at 
first order, through a chiral equation of motion (e.o.m.) which coincides with the usual 
flat expression, dnW = 0. In our situation this cannot be the case as the neglected terms 
come like W, which we suppose naturally large, an therefore no longer negligible. However, 
if the Two-Step procedure is approximately right, the effective theory once the H fields 
are integrated out should be also approximately SUSY. Then, despite the lack of a chiral 
e.o.m. that can be exploited in order to make the procedure fully SUSY manifest, as was 
done in [12], we follow a manifestly SUSY procedure by keeping the auxiliary fields in the 
Lagrangian and integrating out simultaneously all components of the H chiral multiplets. 
This not only leads to a straightforward identification of the effective theory as an approx- 
imate SUGRA theory but also turns out to simplify the integration of the H multiplets, 
at least in a schematic form, in the case where gauge interactions are involved. The fact 
that for the factorizable models the effective theory continues to be a SUGRA one, at least 
approximately, then can be understood from the fact that the weak coupling between the 
two sectors allows to have approximately vanishing solutions for the H auxiliary fields, in- 
dependently of the remaining dynamics. These models present, therefore, an exception to 
the general case studied in [13] where is the hierarchy between the scales what suppresses 
the SUSY breaking effects from the H sector, but we leave a fully superfield understanding 
of this situation for a forthcoming paper [31]. 

The first thing to be noticed, even before turning on the gauge interactions, is that the 
presence of matter fields clashes with the factorizability of the system as their wave func- 
tions in general depend on fields from both sectors and a large mixing can be realized. 
However, if the wave function somehow turns out to be suppressed, the mixing between 
the moduli is still safe, and moreover the dynamics related to the Q multiplets is suppressed 
such that the independence of the H solutions on the other fields continues to hold. In the 



-3- 



LVS this suppression is indeed realized, with the wave functions coming Uke 1/V", with 
the modular weight n a positive number and V the volume of the compact manifold [32]. 
In general one can expect, then, that the corrections to the simplified model be modular 
weight dependent and the reliability of the Two-Step procedure be constrained by this 
number. Nicely enough, due to a necessary tuning on the numerical parameters, and the 
way the matter fields stabilization proceeds, the corrections up to the mass level are not 
only suppressed but also independent of n\. Since for moduli stabilization issues this order 
in the fiuctuations is enough, this universality on the corrections is an appealing feature of 
the LVS as a robust playground for moduli stabilization models. 

Notice that we are being cavalier using the term eff'ective theory: in case there is a mass 
hierarchy we could properly speak about an effective theory in the Wilson sense, but in our 
case things are different and what we call effective theory is simply the one obtained by 
solving the classical e.o.m. and then plugging back the solution on the original Lagrangian, 
resulting with what we will call effective Lagrangian. Nevertheless, the approximate fac- 
torizability of the Kahler invariant function implies also decoupling of the wave functions, 
indeed the scalar manifold is at first order factorizable, so that a two derivative approxi- 
mation is still valid once one requires slow varying solutions in the H sector. This is: the 
fluctuations in the L sector do not excite, at flrst order, fiuctuations in the H sector. 
The introduction of matter fields and gauge interactions in the LVS is far from being just 
and academic exercise. Indeed, the simplest pure moduli realization of LVS turns out to re- 
alize only a deep AdS vacuum, therefore extra ingredients are needed to uplift the vacuum 
being the more natural ones adding extra fields.^ Although developing models realizing 
exponentially sized compact manifold volumina is not as prolific as the ones presenting 
moduli stabilization a la KKLT, mainly due to the rather intricate structure leading the 
stabilization which makes difficult its generalization, there are a couple of proposals leading 
to Minkowski vacua, both of them with matter and gauge interactions [21, 29]. 
The structure of the paper is as follows: section 2 is devoted to generalize the definition 
of nearly factorizable models with the outlook of introducing matter like fields. With this 
generalized setup, but without gauge interactions, we show how the decoupling is realized 
and freezing is a reliable procedure at leading order in a small parameter characterizing 
the mixing in the Kahler invariant function. The analysis is then particularized to the LVS 
where, however, there are in principle several unrelated suppression factors for the mixing, 
namely inverse powers of the volume in the Kahler potential, and non-perturbative effects 
in the superpotential. Besides of recasting and generalizing some of the results of [8, 10] 
the study is also performed in a SUSY manifest way that helps to spot easily the SUSY 
nature of the effective theory and moreover shows itself powerful when dealing with gauge 
interactions; the main target of the paper lies in section 3 where gauge dynamics are intro- 
duced in the study, with the outcome of a constrain on the gauge kinetic function for the 
general setup in order the Two-Step procedure be reliable. The analysis is then performed 
for several LVS models where it is shown that the restriction on the gauge kinetic function 

^Like in the KKLT scenario [17] it is possible to introduce D3-branes at the tip of a throat to do the 
job. This necessarily introduces explicit SUSY breaking terms which one might want to avoid. 
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is naturally avoided. The section closes with few comments about higher order operators 
and its reliability in the simplified model. Also few arguments are given in support of the 
independence of the corrections on the modular weights; section 4 is devoted for the con- 
clusions and some discussion; two appendices are left. One to review the main aspects of 
the LVS models needed for the study, and to introduce a novel model developing LVS. The 
second one gives numerical samples showing explicitly the results obtained analytically in 
the main text. 

2 Only chiral multiplets 

2.1 Nearly factorizable models 

For our purposes it will be convenient to work in a Kahler gauge where the factorizable 
nature of the A/" = 1 SUGRA theory be manifested [8, 9, 30]. This leads us to work directly 
with the generalized Kahler invariant function G = K + Mp log ^^^b-^j where W is the 
superpotential, K the Kahler potential and Mp the reduced Planck mass that for simplicity 
we set to one in the following. The approximate factorizable models are defined as such 
that only suppressed terms mix two sectors identified by H and L, i.e., 

G{H, H, L, L) = Gh{H, H) + Gl(L, L) + eG™,(F, H, L, H) , (2.1) 

with e a small parameter. This form was proposed by Binetruy et al. in [30], to describe 
a SUSY decoupling between the two sectors of fields; however, it lacks from including 
matter like superfields. Indeed, although the form for G can be justified between moduli 
sectors, closed string fields, it is not once matter, open string fields, enter in the game since 
their wave functions depend on both moduli sectors, so factorizability is lost once these 
fields acquire 0(1) VEV. This situation can be easily taken into account generalizing the 
factorizable definition of the system with a form for G inspired on the LVS, where the wave 
function of the matter fields is suppressed by some power of the volume. Thus, in this case, 
splitting the L sector by the moduli, A^, and the matter fields, Q, we propose 

G{H,H ,L,i) = Gh{H,H) + Gm{M,M) + eG.mU<t),4>) . (2-2) 

with (j)^ = {H,Ai, Q}. This is: the matter fields Q only enter in the suppressed part of G 
which can depend in all three kind of fields. This form implies that the matter fields enter 
in the superpotential also in a suppressed way, forbidding 0(1) couplings, e.g., Yukawa 
like. For the moment let us simply say that fields appearing on such a coupling usually 
develop vanishing VEV's, which are not dangerous for the factorizability and we postpone 
its study for the moment. Let us see that a structure like the one defined in eq.(2.2) indeed 
allows for a decoupling between the two sectors by integrating out the H superfields. 
We use the superconformal approach [33-36] to write down the SUGRA Lagrangian since 
the tensor calculus is similar to the one of rigid supersymmetry. The procedure introduces 
new degrees of freedom collected in the conformal gravity multiplet with components the 
graviton, the gravitino and two vector auxiliary fields. Also a compensator chiral multiplet, 
is added. By fixing the scalar and spinor components of the compensator and one vector 
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auxiliary field one recovers the symmetries of ordinary SUGRA. The Lagrangian, then, up 
to two derivatives can be written as an integration over rigid super-coordinates [34-36], 

£ = j de^ (^-3e-^/-'^$^) + j de2$3 ^ ^ ^ (2.3) 

The Berenzin integrals, however, are now deformed by extra-terms with dependencies on 
the components of the gravity multiplet. We will mainly be interested in the scalar compo- 
nent Lagrangian, where these deformations are not relevant, so the calculation is straight- 
forward and easily written like in the global case. With conventions for the metric signature 
(+, — , — , — , — ) and for the chiral fields (p = {(p, Tp, —F^) the scalar Lagrangian reads, ^ 



+ [Gmm - 3 Ga/G,-, J F''^ F'' - WU - Ze' (U + U), (2.4) 

where the gauge fixing <I> = e*^/^, F^ = e'^^^U has been implemented in order to get a 
canonical normalized Einstein-Hilbert action [34] , and we have adopted the convention on 
the notation for the derivatives on the fields BmG = Gm, capital Latin letters running 
over all fields H, Ai and Q. Contrary to the analysis performed in [8, 10] we will keep the 
auxiliary fields in the game. This not only makes the SUSY nature of the theory to be 
completely manifest but also will make the following analysis neat, showing more powerful 
in the case with gauge interactions. Then, to be consistent with SUSY we now integrate 
out simultaneously both components of the H multiplets, whose e.o.m. are 

F':= GiU+(Gif,-^G,Gf^^F^ = 0, (2.5) 
H' := - G.j^d^^'' + G.j^.jj^d^cp^'d^'^^ + GmF^'U + G^j^^F^^'U 

M TpN 



le^/\U + U)G,+ 



GiMN ~ -^{GimG^ + GmG^j^) 



pMpN (2.6) 



with convention of small Latin letters running only over the H fields. The first equation 
together with the on-shell expression for the compensator auxiliary field, U = ^GmF^^ — 
e*^/^, leads to the on-shell expression for the auxiliary fields, that we write as 

Gi = e-^/'G,^F^. (2.7) 

Regarding the approximate factorizability of the function G the e.o.m. for the H fields 
reads, 



G^jF^U + G.jF'U + Gif^F^F'' - Iq^G'^F^f'^ - 1ga/G,jF^^F^ (2.^ 
-^e^/^{U + U)G^ - Gi-d^n^ + Gff^d^H~^df,H'^ = 0{e) , 



^For simplicity of notation, here and throughout the paper, we use the same notation for the chiral 
multiplets and its lowest components, being clear from the context to which one we are referring to. 
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so that at leading order in e the Q fields completely disappear meanwhile the ones 
only appear in the GmF^^ terms. Here we see that although there is no mass hierarchy 
between the two sectors the factorizable nature of the G function, and in particular of the 
Kahler potential, leads to almost decoupled wave functions so that the rotation to normal 
coordinates does not mix the two sectors. In other words, despite the mass scales of both 
sectors are the same the kinetic energy from the L sector cannot, at first approximation, 
excite fiuctuations in the H sector (see also [16]), allowing for a two derivative approxima- 
tion. Looking for slowly varying solutions, using the on-shell expression for the U field and 
eq.(2.5), one gets 

G,,F^U + G.f.F^F' -1{U + 3U)G,-F' - - ^GmG,-F'"f' = 0{e) . (2.9) 

These, have solutions = 0{e) implying, from eq.(2.5), that the lowest component of H 
be solution for 

^^G = 0{e) , (2.10) 
so that once we require that Wh, from Gh = Kh + log | p, be 0(1) necessarily implies 

diGu = ^{diWH + diKHWn) = 0{e) , (2.11) 

which are nothing but the leading F-flatness conditions one finds by solving the e.o.m. 
obtained directly from the lowest component potential [10]. Of course, it might happen 
that these solutions do not exist, or that these do not fix all H fields. We will suppose in 
the following that this is not the case and all H fields are fixed by (2.11) avoiding issues 
like the one pointed out in [37]. The solution for the lowest component of then, can be 
cast as ff* = + eAW, where is the solution for diGu = 0, which is L independent. 
Plugging back these solutions in (2.4) keeping only up to 0{e) factors, one gets, with the 
Greek indices running only over L sector fields, 

^eff = G^-pd^rd^r + G^F^U + 

+ [g^p - \g^G-^ F^F~^ - 3UU - 3e§ iU + U) + 0(^) , (2.12) 

where we have kept explicitly, in the last term, some O(e^) factors in order not to be forced 
to split the G function and overload the notation. Notice that the only place where the 
iJ* lowest component solution appear at C'(e'^) is in the last term. However, expanding 
around the leading solution, diG = 0{e), the corrections will appear only at the ©(e^), so 
we can safely keep only the leading solution for the H^, i.e., HI, as well in the rest of the 
Lagrangian being at most 0{e), and therefore any possible L dependency of the H solution 
completely disappears. 

The effective theory is then described at leading order by the Lagrangian explicitly written 
in eq.(2.12), that is precisely the one of a SUGRA theory with Kahler invariant function 
given by the original one but where the H superfields are frozen at their leading solutions, 

Gsimp{M,M,Q,Q) = G{Ho,Ho,M,M,Q,Q) , (2.13) 
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with Hq the constant chiral superfield with vanishing spinor and auxiliary components, 
and scalar component the leading solution of (2.10). One can further check that the SUSY 
breaking contributions from the H sector are in fact suppressed compared to the L sector 
ones. Taking the canonical normalized F-term, F^^ = \Kj^^j^F^ F^'^l^^"^ no sum, we have 

F^ ~ 0{l) , F^ ~ 0(ei/2) , F'^ ~ 0{e) . (2.14) 

We are finding explicitly that the resulting theory obtained by integrating out the H fields, 
regardless the lack of a hierarchy with the Ai fields, nor with the SUSY breaking scale [13],^ 
can be described at leading order in e by a SUGRA theory whose Lagrangian coincides 
with one of the simplified version. In particular, in absence of matter fields the simplified 
version is valid at next to leading order since for the moduli the Lagrangian is of 0(1). In 
this way we recover the results in [8, 10], with the addendum of making more manifest the 
SUSY nature of the resulting theory. 

We understand, moreover, that it is thanks to the weak coupling between the two sectors 
that the solutions for the auxiliary fields are suppressed, independently of the remaining 
dynamics, and we can realize an effective SUSY theory although the masses for the H fields 
turn our to be of the same order of the SUSY breaking scale. 

2.2 Large volume scenario 

The M = 1 SUGRA theory obtained in type-IIB orientifold compactifications, where the 
LVS are realized, is described by a Kahler potential with general structure, including 0(q;'^) 
corrections [38], given by 

K = -2log{V + ^{S + Sf/'')+ICcs, (2.15) 

where ICcs depends only on the Dilaton, 5, and complex structure, C/*, but whose explicit 
form is irrelevant for our study. The Calabi-Yau 3-fold (CY) volume, V, is a Kahler 
moduli dependent function, such that in case of absence of a' corrections, encoded in the 
^ parameter, the SUGRA theory has a no-scale nature. 

Although it has been studied some generic properties the theory should satisfied in order 
to realize vacua with exponentially large volumina [23], explicit realizations of moduli 
stabilization and characterization of LVS properties are in general based on the so called 
"Swiss-cheese" manifold compactifications which probably encode all the main features of 
the LVS models. In the simplest of such manifolds the compact CY is realized as the 
hypersurface CP^j^ 1 1 6 9]' '^ith /i^'^ = 272 complex structure moduli and /i^'^ = 2 Kahler 
moduli: one breathing mode, T, and a blow-up mode, t. The volume is in this case [39], 

V = A(r^^/2 -t3/2) , (2.16) 

where A = 1/9^2, but we leave it in this form in order to have clearer formulae, and 
we have introduced the notation (pr = (0 + '^)/2 used here after. In the following we 

^Notice that the physical masses, canonicaUy normalized, and SUSY breaking scale, once the L sector 
breaks it, are all of C(l). 
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stick to this form for the volume, mainly to match the existent models but also to have a 
precise dependency on the L fields being crucial for the analysis. Our results, however, are 
automatically extended to more complicated Swiss-cheese manifolds since in general with 
more blow-up modes the volume is generalized as V ~ J'S/^ _ ^ . t^J"^ . At the qualitatively 
level our results are also expected to hold for other kind of manifolds. 
There is a flux induced tree level superportential given by [40] 

W,s = Jn^AGs, (2.17) 

where is the (3,0)-form of the compact manifold, and ^3 = ^3-!- iS H3 the combined 
three-form-flux. The Kahler moduli, instead, appear in the superpotential only through 
non-perturbative dynamics, which can also depend on the other moduli. In the spirit of a 
Two-Step moduli stabilization the Dilaton and Complex structure moduli are regarded as 
frozen due the superpotential (2.17) with larger dynamics compared to the one encoded in 
the non-perturbative part. It was found in [10, 12] that in order this procedure be reliable 
for generic Kahler potential one needs to require a tuning at the vacuum for the tree level 
superpotential. Indeed, this tuning implies a mass hierarchy which warranties the decou- 
pling between the two sectors. In a more natural scenario the classical superpotential is 
to get 0(1) VEV, and the decoupling relies in the factorizability of the G function [8, 10]. 
On the side of the superpotential the mixing between the complex and the Kahler moduli 
appears only by non-perturbative dynamics, so are naturally suppressed and the approxi- 
mate factorizability, as defined in eq.(2.1), is safe. On the other hand, the a' corrections of 
the Kahler potential, breaking the factorizability, are not naturally suppressed. However, 
whenever we stand at points of the moduli space such that the compact manifold volume is 
quite large, like in the LVS, these corrections are small an approximate factorizable system 
is realized. Notice that contrary to the KKLT like models studied in [10, 12], besides the 
suppression factor coming from the non perturbative dynamics, there is a second one dic- 
tated by the size of the volume. This two a priori unrelated factors are however connected 
through the stabilization of the moduli and this relation depends on the particular model, 
making not viable a generic study as was possible for the KKLT like family of models. 
As said before the presence of matter fields breaks the factorizability in the Kahler po- 
tential, but in the extended sense given in eq.(2.2) one can still get near factorizability 
between the H sector. Complex structure and Dilaton moduli, and the L sector, Kahler 
moduli and charged fields. Indeed the Kahler potential for the matter fields reads [32] 

KD^\Q\\ (2.18) 

where Q is a generic charged field, omitting for simplicity possible index contractions, Z is a 
real function of the Complex structure and the blow-up Kahler moduli t, and the modular 
weight n is a positive number. Depending on the cycle the D7-branes wrap, where the 
matter is realized, the moduli dependent function Z trivializes or not, and the suppression 
factor is a power of the volume or just a power of the breathing mode T. For our purposes 
this schematic form is enough, noticing again that around points in the moduli space where 
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the size of the volume is quite large the appearance of the Q in the generalized function 
G is suppressed as far its dependencies in the superpotential be also suppressed. We have, 
then, in general a third suppressing factor depending on the modular weight. 
Let us start by considering the simplest case where W \s Q independent, so that the 
superpotential takes the form, with A a function in general of the H fields, 

W = Wsc + Ae-''^ (2.19) 

with the non-perturbative part of the superpotential depending only on the a priori re- 
garded small modulus. More precisely the breathing mode is supposed extremely large so 
that such dependencies vanish. If we now suppose that the Kahler moduli stabilization 
follows like in the pure moduli case, as actually does and even it turns out to happen in 
realistic models seed by this one (see appendix A), the non perturbative suppression factor 
coincides with the inverse of the volume and we can define a single suppressing factor given 
by e ~ 1/V ~ l/r3/2 ~ Ae""*, meanwhile t = 0(1). The analog of eq.(2.9) then reads 

G,,rU + G,f,rF' - ^G,,Gj^F^F^ - ^GijGMF^F^ '1^^ + ^U)GfjF~^ 

^{U + F^ + G^F^)Gi,F- + G.^^F'^F^ , (2.20) 

where as before the capital letters, M and N run over all fields, i, j, over the H ones, and 
the Greek indexes over the Kahler moduli and the Q fields. In the r.h.s. we have been 
sloppy with the holomorphic indexes and their contraction being only interested in the 
leading scaling in e: for example for t, in case Q has a vanishing VEV, we have G^^ ~ 
meanwhile Gtn ~ so we take the last one for the analysis. The vector and matrix shown 
in the r.h.s. of eq.(2.20) have the scalings on e, 

G'.«~(e2/3(e + e"Q2)^^^^ng2^^ng 

/e4/3(e + e"Q2) ,2/3(^2 ^ gng2^ 

e2/3(e2 + eng2) ^ + e^\Q\^ e"Q | , (2.21) 

^ g2/3+ng ^ng 




i fixed and choosing the ordering {T,t,Q}. Having in mind that the on-shell solution for 
the auxiliary fields in the L sector are necessary suppressed by powers of e and looking for 
approximate SUSY solution in the H sector the leading solution for F^ auxiliary fields has 
the schematic form 

F' ~ Gi„F" ~ £2/3 ^ er^g2)^T ^ ^ gng2^ pt ^ ^UQpQ ^ (2.22) 

being again sloppy with the holomorphicity of the indexes, in fact the solution will mix 
them up in general. This leads us again to conclude that the leading solution for the H 
lowest component is dictated by the the F*-fiatness conditions, 

e^/'^Gi = Gif^F^ ~ e2/3(e + e"Q2)^T ^ ^ e«Q2^^i ^ ^uQpQ ^ ^2.23) 

then, although e*^/2 ^ ^j^e on-shell expressions for the show that Gj vanishes at 
leading order in e. Indeed equation (2.22) tells us about the suppression of the SUSY 
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breaking contribution from the H sector. As before, we consider that Gh = /Ccs+log |VFscP 
is such that ahows for such solutions and fixes the whole set of H fields. Plugging back 
these in the Lagrangian one has 

= G^^d^rOf^r + GaF'^U + GaF^U 

+ (g^-^ - lOaG^^ F^f'P - 3UU - 3e^{U + U) + 0{{r)^) , (2.24) 

where as before we can keep only the leading solution for the scalar component of H, i.e., 
diGu = 0, which are independent of the L fields, being the further corrections negligible. 
As pointed out before the corrections 0((-F*)^) will have not only the SUSY combinations 
FF but also some pure holomorphic and pure non-holomorphic terms that reveal the non- 
SUSY nature of this effective Lagrangian. We should compare all the 0((-F*)^) corrections 
with terms coming from a SUSY Lagrangian constructed from Gsim{T,T,t,t,Q,Q) = 
G{HQ,HQ,T,T,t,t, ,Q,Q), with Hq a superfield with vanishing spinor and auxiliary field 
components and whose lowest component is fixed by the leading solution to the F*-flatness 
conditions, which are again L independent. This Lagrangian is precisely the one explicitly 
written in (2.24) with some subleading corrections due to the fact that the correct solution 
for the scalar component in shifted from the leading SUSY value. 

Let us perform explicitly the analysis of the term |-F*P, so to make the procedure clear, 
and furthermore spot possible numerical factors necessary to precisely understand the 
numerical tests performed later on in the appendices. Taking the term we have that 
in the simplified version it is lead by 



This term is also generated in the effective Lagrangian by \GitF^\'^, with leading terms 
Git ~ o e -^'d^A+^\Q\^ ~ e+e"|gP, so there is a further contribution not present in the 
naive Lagrangian which goes like (e + e'^|Qp)^|F*p, but is suppressed by A = 0{e, e"|(5p). 
In order to keep track of the numerical factors one has to recall that Ae^^^ = AX\/i ^^^°^ e 
(see appendix A), with A encoding possible extra factors missed by the analytic procedure, 
so that in the Q independent parts of the corrections one finds an enhancement. Being 
precise one finds A ~ f At3/2i2g^ Tj^g 

same analysis is done for all terms, which in the 
simplified Lagrangian are ruled by the derivatives 



G'a~(e2/3(l + e"|Q|2),e + e"|Q|2,6-Q 

/e4/3(l + e"|Q|2) e2/3(g + gn|g|2) g2/3+ng\ 

e2/3(e + e«|Q|2) e + e^\Q\^ e"Q 



(2.26) 



In order to compare terms that are not originally present in the simplified version, e.g., the 
pure holomorphic {F'^)^ term, we replace one of the auxiliary field to its scaling around 
the vacuum, which from eq.(2.26) with ~ e'^/^G~^Ga no sum and (2.22), are 

F^~ei/3, F*~e, ~ eQ , ^ F^ ~ + e'^+iQ^^ (2.27) 
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so that 



[F'f D (e^ + e"+iQ2^ (^e'^/^e + e^Q^)^^ + (e + e''Q^)F' + e^QF'^) , (2.28) 

and compare the results with the corresponding one, e.g., i^^, regarding the scaling for 
the compensator auxiliary field as C/ ~ e, resulting from its on-shell expression, U = 
\GmF^ -e^l'^. In this way one finds that all corrections encoded in the 0((-F*)^) terms 
are suppressed by 0(e,e"(5^). As explained in appendix (A.I) in order the stabilization 
of the moduli indeed proceeds as in the standard case and we can relay on the given 
scalings, one should impose the condition (Q^) < e^~", so the largest corrections are of 
0{e), independent of the modular weight n!. In this way we generalize the results obtained 
in [10] working now in presence of matter fields and keeping track of the auxiliary fields. 
To close this section let us give a look at the canonical normalized auxiliary fields VEV, 
showing that indeed the SUSY breaking contribution from the H sector is suppressed, 

Fj~6, Fi~e3/2^ ~ ei+"/2Q < e3/2 ^ F^, ^ e\ (2.29) 
as obtained from eqs.(2.22), (2.26) and (2.27). 
3 Chiral and vector multiplets 

We now turn on gauge interactions, that can be seen as some gauged isometries of the scalar 
manifold, which for a A/" = 1 SUGRA is of the Kahler type. The chiral multiplet trans- 
formations are then defined by 6(f)^'^ = K^X^ , 5cj)^^ = A^A^^, with A^ the gauge chiral 
parameter, (X^), the holomorphic (antiholomorphic) Killing vectors generating the 
(gauged) isometry group Q, and the gauge group indices running over A = 1, 2 . . . adj(t/). 
The associated gauge vector fields transform as dV"^ = —i{A'^ — A ), inducing a total 
transformation on G of the form 6G = A^A^G/ + A'^A^Gj - i(A^ - A^)Ga = 0, with 
Ga denoting derivatives of G with respect to the vector multiplets, and as before the Latin 
letters running over the chiral multiplets. Gauge invariance of the system then reads^ 

GA = -iXiGi, (3.1) 

a relation that we will use in the following. Notice that if we want to fix the H fields 
by the conditions Gi = 0, all H fields should be necessarily neutral otherwise by gauge 
invariance, eq.(3.1), the equations turn out to be linearly dependent and leave some unfixed 
directions.^ This requirement also avoids the possibility the H fields be sourced by the 

*We have not included constant Fayet-Iliopoulos terms in eq.(3.1) since they seem not to appear in 
theories raised from string compactifications. 

^Notice that one might generahze the way the H fields are fixed by regarding two gauge sectors, say 
Q — (g) C/A, so that the H are charged under Q ^ but neutral under Qa, and the opposite for the L 
fields. Then fix the H fields by the requirement Gi = and = 0. This is, besides the F-flatness the 
/3-flatness conditions. In this case some of the H field are then stabilized by gauge dynamics and not by 
the superpotential, as usually regarded for the H sector in string compactification, so we left this situation 
out of our analysis. 
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gauge fields in the effective theory. We set, then, in the following = 0. 
Using the conformal formalism the Lagrangian of eq.(2.3) is first of all corrected by the G 
dependency on the vector multiplets V"^, promoting to covariant derivatives the derivatives 
in the kinetic terms, i.e., 

9^0' ^ d^(t>' + XiV^ . (3.2) 

Further terms for (2.3) come from the chiral Lagrangian and from the kinetic part of the 
vector multiplets Lagrangian. Working in the Wess-Zumino gauge, denoting the vector 
auxiliary fields by one has the following extra terms to the ones in eq.(2.4) [33] 

C D GaD^ + hiABD^D'' , (3.3) 
with hAB = ReifAs), Jab the field dependent gauge kinetic functions. 
3.1 Nearly factorizable models 

As before let us first tackle the case of the generic nearly factorizable models where the 
situation follows as before with the e.o.m. for auxiliary fields not changed, eq.(2.5), and 
for the scalar components the analogous of eq.(2.8) now reads, 

-Gfjd'^W + G^jid^Wd^'H^ + GijFW + G^f^F^F^ - ]^{U + 3U)GijT 

-^GijG^F^F'' - ^GmG.jF^'T + ^^D^D^'d.hAB = 0{e) ,(3.4) 

obtained by using the relation (3.1), the neutrality of the H fields, so that OiGa = 
—iX'^Gai, and regarding that none of the Killing vectors X" get anomalous large val- 
ues, i.e., O(X^) ~ 1. We have, then, that the presence of the last terms makes the solution 
~ e no longer automatic; indeed, the e.o.m. tells us that whenever the D-term gets a 
non-vanishing VEV, which is naturally of 0(1), it back reacts on the H sector, even if neu- 
tral, inducing a large F*-term. In particular the auxiliary fields solutions would depend 
strongly on the L fields, therefore cannot be simply neglected, as a Two-Step stabilization 
requires, and have to be properly integrated out. Notice, moreover, that the mixing be- 
tween the F and the D auxiliary fields makes again explicit the non-SUSY nature of the 
effective theory, which on the scalar side still is a two derivative description.^ 
The decoupling would be still safe if somehow the dependency of the gauge kinetic function 
on the H superfields is absent or suppressed by say e. In this case things turn out to be like 
in the case without gauge dynamics as the solution ~ e is again valid, and the leading 
solution for the lowest component is again dictated by the leading F-flatness conditions, 
that do not depend on the L sector. Thus, the simplified theory, where the H superfields 
are regarded as frozen at their leading solution, is reliable at leading order, as far as the H 
superfields are neutral and the frozen value is, up to 0{e) corrections, the leading solution 
to the F*-fiatness condition. 

Notice that we have not distinguished between the cases with and without breaking of the 
gauge symmetries. In the study done in [12] it was important to keep track of the gauge 

®A similar analysis was done in [13] studying the SUSY description of effective theories. 
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symmetry breaking as some of the fields in the L sector, regarded as hght, acquire heavy 
masses via the Z)-term dynamics. Then, one needs to distinguish these fields from the 
low energy spectrum when speaking about an effective theory integrating these fields out, 
not being in general possible to freeze them out, together with the H fields. In our case 
however, modes acquiring masses through this mechanism are not heavier than other fields 
in the L sector. Indeed, regarding the eigenvalues of Hab of 0(1), their masses should be of 
the order of the eigenvalues of vector field mass matrix associated to the broken generators, 

mIb = 2GjjX'aXb. (3.5) 

which are clearly of 0(1). Therefore, there is not need of spotting such fields and no formal 
distinction appears if the symmetry is broken or not. There might appear a small hierarchy 
with such fields, like in the case Nf < Nc worked out in section A. 2. In this case one can 
proceed like done in [12] by integrating the massive vector superfields through the superfield 
vector equation dyK = [13] and choosing a convenient gauge fixing, fixing the value of one 
of the charged fields with a non-vanishing component in the would-be Goldstone direction. 
Since the Kahler potential has a factorizable form and the H multiplets are neutral the 
dependency of the solution on the H fields will be necessarily suppressed. Plugging back 
the solution into the original Kahler potential leads then to a new Kahler potential which 
however has still a factorizable form, and the same analysis can be performed with only 
the unbroken symmetries. However, since in our models is the weak coupling between the 
sectors and not the mass hierarchy what plays the important role in the decoupling, these 
issues are not relevant for our study, being this more clear if we remember that the effective 
theory we are speaking about is not a proper one in the Wilson sense. 
We find, then, that in order Two-Step procedure be reliable in general^ for this generic 
factorizable models, a further restriction is required, namely a suppressed dependency of 
the gauge kinetic function on the fields to be frozen, something that cannot be justified in 
general for realistic situations. Nevertheless, it is clear that this general setup as stands so 
far cannot be connected to a real situation as SUSY is broken at Planck scale. Interestingly 
enough, in the next section we find that in realistic realizations of nearly factorizable 
models, namely the LVS, this restriction is avoided naturally. 

3.2 Large volume scenarios 

In type-IIB orientifold compactifications with Z)7-branes the role of the gauge kinetic func- 
tions is played by the Kahler moduli parameterizing the size of the cycles where the branes 
wrap. The presence of fiuxes, however, induces corrections on the gauge kinetic functions 
that depend on the Dilaton modulus [41], which are supposed to be frozen in the Two- 
step stabilization spirit. Thus, the condition we found for a reliable decoupling in generic 
nearly factorizable models is not naturally realized. Since the mixing between the and 
the D"^ auxiliary fields, which is what rises this problem, still appears, if decoupling indeed 
applies its justification should come from other means. In the previous case, however, we 

^We will show in later sections that what really matters is the VEV of the , so decoupling is still 
save whenever this VEV is tiny. 
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were forced to strengthen the requirements on the models from the fact that the are 
naturahy of 0(1), as seen by their on-shell expression, 

j^A ^ _hABQ^ ^ ih^^X'XGo, , (3.6) 

with the matrix h"^^ = h'^^, since apart from the matter field contribution, of 0{e), 
D"^ Ga = 0(1). Even more, this scaling naturally holds at the vacuum, where the 
following dynamical relation is satisfied [42], 

J (m1^ + hAB{GjjF'FJ - e^)) = G^jjF^F'^, (3.7) 



2 

everything being naturally of 0(1). On the other hand, for our type-IIB setup all G^ are 
suppressed by powers of the volume, so in a LVS things seems to be controlled on favor of 
the decoupling, but even though the decoupling is still not trivially clear. Moreover like in 
the pure chiral fields case a complete generic analysis is not possible, so we will follow the 
study separately for explicit models with matter fields realizing moduli stabilization in LVS 
with vanishing cosmological constant [21, 29]. In this models the main effects from gauge 
dynamics, besides the non-perturbative superpotential, come from Abelian sectors under 
which the moduli get charged, so in the following we only consider tree level effects coming 
from such U{l)x group. This do not loose the generality of the study as the analysis 
of other sectors follows in the same way and even more straightforward as the moduli are 
neutral under those sectors and the dynamics are simpler. As a matter of fact it is common 
to neglect such non- Abelian sectors by working in a meson field description following D-Hat 
directions. 

As usual the symmetry is linearly realized for the matter like fields, instead the moduli, 
controlling the gauge kinetic functions, can develop a non- linear realization due to a Green- 
Schwarz (GS) mechanism that cancels pseudoanomalies of the Abelian sectors. In the case 
of type-IIB superstring with D7-branes the Kahler moduli get charged once world volume 
magnetic fluxes are turned on.® The U{l)x symmetries in the chiral superfields are then 
described by the holomorphic Killing vectors, 

= i<lx<t^^ no sum , Xx = i^Sx' , (3.8) 

with Qx the charge of the chiral superfield (f)^ and the GS coefficient a real parameter 
associated to the modulus . The gauge invariant Kahler potential for the moduli has 
then the functional dependency 

1 

— ( 
2 

with the vector superfield associated with the U{l)x symmetry. To make simpler the 
analysis we will regard that only one of the Kahler modulus gets charged, but in general 
is possible that several of them develop non trivial charges. The main effect for moduli 
stabilization issues is the new contribution induced in the D-teim dynamics which enters 
like a Fayet-Iliopoulos (FI) term but now field dependent [21, 29, 45-49]. 



(3.9) 



*To be precise a topological condition should be realized, namely that the 4-cycle whose volume is 
characterized by the Kahler modulus has an intersection with the 2-cycle where the world volume flux is 
non-trivial, in order such a charge be not null [43, 44]. 
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3.2.1 Charged large modulus 

The simplest implementations of this kind of models were developed in [21] with the large 
modulus charged.^ Although the explicit realization of the models might have several 
matter fields for our purposes their final out-shot can be cast using only one Q, so that 

X],=i^6^, X*,=0, X^=qQQ. (3.10) 

By gauge invariance the superpotential does not depend on the Q field so the F-term part 
of the scalar Lagrangian goes like in the case without gauge interactions, and it is possible 
to stabilize the Kahler moduli like in the pure moduli case. The scalings shown in eq.(2.26) 
and eq.(2.27) then hold, and new scalings needed for the analysis are, using eq.(3.1), 

Gx ^ hx D"" 6Te^/^ + e^'lQW Gxi = -iXf ~ ^te'/' + • (3.11) 

Thus a small size for is implied and again is possible to find a suppressed solution for 
the auxiliary fields, using the on-shell expression for the auxiliary fields, 

F' ~ + -GxD^ 

e 

= e2/3(e + e"Q2)^T ^ ^ e"Q2^F* + e"QF« + (^Te"^/^ + e""i|Q|2)D-'^ , (3.12) 

with the last term 1/e factor coming by regarding the compensator scaling [/ ~ e. This last 
term apparently couple strongly the two sector, and might give rise to L dependencies in 
the H solution; however, it turns out to be inoffensive as explained bellow, so the effective 
Lagrangian is again given by the simplified version plus corrections of the form 

OiF'f ~ (F>)2 + -GxD^Fj, + \g\{D^)^ , (3.13) 

where by Fp we denoted the pure F auxiliary fields contribution, eq.(2.22). 
The analysis of the corrections to the simplified version from the {Fp)"^ part is exactly the 
same as the one with no gauge interaction since the F-term part of the Lagrangian does 
not change. Remains, then, only to check the new features coming from the last two terms 
which enter as corrections to the F", and {D-^Y terms. 

Under this setup there are two possible scenarios, depending on the sign of the Q charge, 
differing by the VEV of Q which is ruled mainly by the minimization of the D-term 
dynamics, eq.(3.3), which once the D auxiliary fields are integrated out, eq.(3.6), reads 

where the gauge invariance of the superpotential, X^djW = is used. 

The leading behaviour of the first term in (3.14) is given by OtK ~ —l/T^ < 0, and if 

the charge is of the opposite sign of the GS coefficient the minimization of Vd leads to 



^Although the authors of [21] also study models with a charged small modulus, these turn out to lack 
of realizing a vanishing cosmological constant. 
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a vanishing VEV for Q. Then, at the minimum the D-term potential scales like Vd ~ 
7^ (f^) ' ^^"^ ^° order to uplift the Vp = —0{e^) vacuum without sweeping away 
the moduli stabilization, one requires a tuning 5'^ ~ e^^^^Jhx- This, together with the 
vanishing of (Q), implies that the dangerous coupling in (3.12) is in fact suppressed even 
in the extreme case when hx ~ T,. ~ e"^/"^. From the effective Lagrangian point of view 
we have that since GxD^ /e ~ e^, independently of the scaling of hx, and Fp ~ as 
before, then GxD^ Fp/e leads to correction to the pure F-term part similar to the ones 
coming from (Fp)'^ and which we learned before are suppressed by 0{e). These also rise 
corrections to the GxF>^ term of the form eOxD^ . The term hx{D^)^ is corrected by 
the last term in (3.13) that using (3.11) and the constrain on (5^ is of the form ehx{D^)^ . 
We find, then, that all new possible corrections continue to be suppressed at least by 0{e). 
In case the charge turn out to be of the same sign of the GS coefficient the minimization of 
the D-term potential tends to a cancellation lead by a non- vanishing VEV for Q, which now 
scales like (Q^) ~ e^/s-"^^. The uphfting now proceeds from the F-term part with a term 
scaling like Vp D e^"''"Q^, so in order to uplift without spoiling the moduli stabilization 
a tuning 5^ ~ eV3 is required. Due to the Q dependency in the F-term potential the 
cancellation of the D-tevui part is not exact, finding in fact that at the vacuum ~ e^, 
independently of the scaling of hx, and therefore Gx ~ hxe^- We have then that GxD^ ~ 
hxe^ and even in the worst case hx ~ T ~ e~^/^ it is smaller that in the previous case so the 
correction to the pure F-term part are further suppressed, more precisely are suppressed 
at least by ©(e''/^), the same happening to the corrections to the hxiD^)"^ terms. Since 
the change in Fp is negligible, Fp ~ + hxC" ~ e^, the corrections to GxF>^ are still 
of 0{e).^^ Interestingly we still find that the corrections are independent of the modular 
weight and the decoupling is in general safe. 

3.2.2 Charged small modulus 

Once one allows a charged small modulus things turn out to be more reach and intricate as 
gauge invariance imply the appearance of matter like fields in the superpotential. Thus the 
stabilization of the moduli should be revisted, possibly affecting the scalings used above. 
Again the minimization of the D-term potential, which now takes the form 

Vd-'j^ [-^dtK + QQQdQK^ ~ (-5*6 + gge^igp) , (3.15) 

plays a crucial role in the stabilization of the matter. Notice, however, the change in 
the sign for the field dependent FI term which makes the implementation of such setup 
quite different to the usual SUGRA FI model [45-49]. A brief review of models with these 
characteristics is done in appendix A. 

^°One might wonder if it is correct to take the scaling of at the vacuum, which for this last case is 
drastically smaller than the off-shell one. In fact, taking the off-shell scaling (3.11) with the constrain on 
one would find corrections of 0{e^^^) /hx so that in the situation hx ~ f ~ 1 the corrections would be 
C'(e^''^) not of 0{e) as claimed above. We have check numerically that indeed this is not the case by taking 
a sample with hx = t,., finding still corrections of C(e). 
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Np < Nc case 

In the window Np < Nc, with Nf the number of famihes and Nq the rank of the gauge 
group, that we suppose to be an SU{Nc), the system develops an ADS superpotential [50] 
which in the particular case Np = 1 and Nc = 2 takes the following form 

„— at 



W = Wcs + A- 



(3.16) 



where the chiral matter fields charged under the non-Abelian sector are described by the 
canonical normalized meson 6 - 



2QQ. As before Wcs and A are 0(1) functions of 
H fields. We have also added a mass term, in general i?-dependent, which from gauge 
invariance under the Abelian sector depends on a second matter like field p. We normalize 
the charges such that (5* = 2/a, thus = —1/2 and qp = 1. This model although proposed 
in [21] was not studied there, so we take the opportunity to carefully check if indeed leads 
to a LVS and characterize its vacua (see appendix A). Despite the fact it turns out that 
as it stands it cannot lead to realistic solutions with vanishing cosmological constant we 
do the analysis of the decoupling having the model interesting characteristics like the mass 
term in the superpotential. The modular weights for both matter fields are taken equal and 
the kinetic gauge function is taken to be fx ~ t- The VEV of </> is ruled by the cancellation 
of the D-term dynamics, meanwhile for p a non- vanishing VEV is only possible for a non- 
zero mass parameter. The moduli are ruled by a potential whose structure is similar to the 
standard LVS but with different powers of the volume. Nicely enough a fully analytical 
solution is found in the large volume expansion, with out-shot. 



-at ^2- 

e ~ e 



t ~ 1 



e 2 



(3.17) 



Contrary to the first two cases, here it is not necessary to argue for a tuning in the GS 
coefficient, however, the mass term may lead the dynamics and sweep the solution, so 

3n— 1 

should be tuned by m ~ e^~. With this consideration the following scalings are found. 
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with ordering {T, t, (p, p}, and where for simplicity in the analysis, contrary to the previous 
examples, we have replaced the matter fields by its VEV to extract the scalings. The mixed 
derivatives scale as 
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from which we estimate the VEV of the auxiliary fields, 



~ £i/3 ^ e, F'^ ~ , ~ and F' ~ e^. 



(3.20) 
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We have also Gx ~ ~ e^. The canonical normalized auxiliary fields VEV's, encoding 
the SUSY breaking contribution from each sector, scale as FJ ~ e, 

Fc ~ e'^/^ and ~ e^. The same analysis followed before shows corrections to the simpli- 
fied version suppressed by 0{e^^^). Nicely enough, even though we are now dealing with a 
fourth suppression factor, namely the mass parameter, the leading order for the corrections 
to the simplified version are again independent of the modular weight. Like in the other 
cases, this is due to the tuning on the parameters from the requirement of a vanishing cos- 
mological constant, or more weakly, the requirement that the stabilization of the moduli 
is not spoiled. A related result is the independence on the modular weight of the VEV 
of non-pertubative superpotential, being in all cases ©(e), and even the mass term in this 
last example goes like mpcfp' ~ e^/^. 

The largest corrections we are spotting here come from Gip ~ (fp'dim ~ together 
with GiT and Gu affecting the GxpF"^ F^ + GtpF^F^ + h.c. terms in the simplified version. 
Therefore, in case the mass parameter is absent, or independent of the H fields, the correc- 
tions are smaller and turn out to be like in the previous cases 0{e). In the numerical check 
done in appendix B, for example, it turns out that the 0{e) corrections are as important as 
the ones coming from the mass parameter due to enhancement factors like the ones found 
before. 

Nc <Nf < 3/2Nf case 

A realistic model already present in the literature is the Krippendorf-Quevedo model 
developed in [29], regarding the Seiberg dual description of the theory in the window 
Nq < Np < 3/2Np. The model can be minimally described by two matter fields with the 
following superpotential, 

W = Wcs + Ae-"^ pel) , (3.21) 

and as before Wcs and A can depend on the H fields and are reg arded to be 0(1). The 
charges can be normalized fixing the GS coefficient to be (5* = 2/a so that, = — 1 and 
Qp = 2. The same modular modular weight is taken for both matter fields and the kinetic 
gauge function is taken to be /x ~ i- Like in the previous case the VEV of (j) is ruled by 
the D-term minimization. A rough analytical approach shows the scalings to be 

e-«*^e"+i/2, tr^l, 0~e{i-«)/2, p ^ e{2-«)/2 . (3.22) 
Then for the analysis we have the scalings 
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^^See appendix A. 2 where is argued that such assumption breaks down whenever the volume turns out 
to be huge. Here we simply disregard such subtlety which makes obscure the analysis. 
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with ordering {T, t, (j), p} and where hke in the previous case we replaced exphcitly the VEV 



of the matter fields so to make simpler the analysis. In the vacuum also D-^ 
which together with the mixed derivatives with the H sector 
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(3.24) 



imply 



4 — n 
~ g 2 



(3.25) 



The analysis follows like in the previous cases finding corrections suppressed at least by 
0{e), lead by the equivalent ones found in the pure moduli case; indeed, the VEV for the 
non-perturbative part of the superpotential again scales like 0{e) and the corrections turn 
out to be independent of the modular weight. Again the canonical normalized F-terms: 



3/2 ^ ^3/2 pP 



~ e 



and Fl ~ 



show the suppressed contribution 



from the H sector to the breaking of SUSY. Unfortunately the rough analytical approach 
we followed in finding the solutions is not enough for our computer skills so to allow us to 
find numerically this kind of vacua and check the conclusions just raised. 

Higher order operators and universality in the corrections 

So far we have regarded only the lowest order operators in the Kahler potential neglecting 
contributions of the form K D Q^/V" with p > 2. Although the introduction of such terms 
implies to deal with extra suppression factors, as the power m is in principle unrelated to 
the modular weight n, we can from our results easily spot the effects of integrating out the 
H fields in such a case. In section 2.2 we found that the simplified version misses terms 
that are suppressed by 0{e^Q'^), where the can be seen as fluctuations of the matter 
fields. This means that the full effective theory realizes terms like Ceff D e^'^'^^Q'^, which 
affect terms coming from the higher coupling in the Kahler potential with p = 4, namely 



■'Simp 



2+m 



IQI^, e.g., from Gqq\F'^\'^ , telling us that such a term in the Kahler potential 
is not reliable unless m < 2n. Similar comments follow for any other higher couplings. 
These higher order effective terms can be easily understood by collapsing Feynman dia- 
grams with H fields in the internal lines, finding also the same kind of corrections once 
higher order operators are regarded in the superpotential; then, a 0(1) Yukawa coupling 
can induce operators correcting the operator in the superpotential and the one in 
the Kahler potential, exactly as is found even in case a mass hierarchy is realized [12]. 
The above discussion is somehow related to the fact that contrary to the moduli, the 
derivatives with respect the matter fields always carry a decreasing on the suppression fac- 
tor controlling the factorization of the theory, since these fields scale with inverse powers 
of the volume. Then, at some order in the matter fields fluctuations the factorizability is 
expected to breakdown and the corrections start to be non-suppressed. 
Notice that this analysis points out also, and contrary to claim done above, for the possi- 
bility the corrections to the lower order operators to be dependent on the modular weights. 
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Indeed, the corrections of 0{e"'Q'^) can be seen also as 0{e^Q)6Q corrections with 6Q the 
fluctuations. Let us see more precisely that in fact up to the mass level, at least for the 
models studied, the corrections are universal. 

For this we write down the scalar potential expanding in fluctuations 6Q and 5H around 
the Ho solution, 

V = Vo + V^5H + Vq5Q + VQi6Q5H + Mfj5H5H + Vqq6Q5Q 

~ + e^6H + e^JQ + e^6Q6H + e^6H6H + e^+'^dQSQ , (3.26) 

where we have used the scalings at the vacuum, and the fact that Vqi ~ Vq ~ QVqq, valid 
for the fields not stabilized by the D-term cancellation in which case a larger suppression 
appear in the Vq and therefore things are even safer. We also used that the non-canonical 
mass squared for these matter fields turn out to be always of 0{e'^^^) and the largest 

1 — 71 

VEV for the matter is 0{e'^). An expansion in the moduli fluctuation is also implicit in 
the coefficients but for our present purposes are not relevant. Then an integration at the 

l + n 

gaussian level leads to 5H ~ e + e 2 SQ, and plugging these back in the potential we have 
corrections of the form 

5V = e^ + e^5Q + e^+''6Q5Q, (3.27) 

which are suppressed compared to the original operators by 0(e), independent of n. These 
as far the largest corrections are concerned, since we took the largest VEV for the Q's. 
This analysis, in fact, is quite general as uses the largest possible VEV for a generic Q field, 
in an arbitrary LVS model, therefore the conclusion is expected to hold in a more general 
context. 

4 Conclusions 

Our results show how the factorizable models present a framework where freezing of a set 
of light fields is a reliable procedure. Although several papers have worked on these models 
[6, 8, 9, 14], we present a generalization which describes the introduction of matter like 
fields, whose wave function depends in general on fields from both sectors. Then, once 
these fields enter only in the suppressed part of the generalized Kahler function, eq.(2.2), 
decoupling between the two sectors is easily understood through the same analysis done in 
the pure moduli case, with the same conclusion: freezing of the H fields is reliable as far 
as the H lowest components stand at their leading F-flatness solutions. The corrections, 
however, contrary to the pure moduli case, where start at the next-to-next-to-leading 
order in the suppression parameter e, start now at the next-to- leading order. With gauge 
interactions turned on this conclusion does not change but one has to require the H fields 
to be neutral and, in case the D-term SUSY breaking be non suppressed, that the gauge 
kinetic function dependency on these fields be suppressed. These models although not 
realistic, with a SUSY breaking scale or order the Planck mass, present all the important 
features needed to understand how the freezing of the Dilaton and Complex structure 
moduli proceeds in the LVS of type-IIB string compactifications. Nicely enough, in these 
scenarios the D-term SUSY breaking is naturally suppressed by some powers of the volume 
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and the constrain found for the gauge kinetic function is avoided. 

A somehow unexpected result is that, in ah models explicitly studied, the largest corrections 
to the simplified version, up to the mass level in the matter fluctuations, are independent 
of the modular weight for the matter fields. This due to the tuning on the parameters and 
a constrain on the largest possible VEV a Q field can get, both related to the requirement 
of a vanishing cosmological constant, or the more relaxed and general one of avoiding the 
destabilization of the vacuum. On the other hand, higher order corrections may depend on 
the modular weights. These, however, are expected to be irrelevant for moduli stabilization 
issues and therefore harmless for the reliability of the Two-Step procedure. 
Our study for the LVS, although far from being fully generic, shows how the Two-Step 
stabilization procedure is expected to be safe in realistic and complicated scenarios, and 
how the corrections are universally independent of the modular weights. It also introduces 
a simple way to estimate the mistake done by performing quantitative analysis from the 
simplified version. Still, a more conservative approach can use the fact that in moduli 
stabilization studies the harder part for the quantitative information is finding the vacua, 
say numerically. Then, once the vacuum is found through the simplified version one can 
implement the full model knowing that the solutions are just slightly perturbed. This, 
in fact, is the procedure followed in the analytical and numerical examples shown in the 
appendices. With the full information at hand one is safe of missing fine information 
probably needed for precision tests, for example soft terms ruled by the H fields. Off course, 
explicit string implementations most likely would need a sort of Two-Step procedure, even 
if less drastic, due to the huge number of fields. Nevertheless, our study sets down a 
consistency check for the LVS as a systematic procedure for moduli stabilization. 
An interesting result our analysis shows, is the SUSY nature of the effective theory despite 
the lack of a hierarchy between the energy scales in both sectors, in particular with the 
SUSY breaking scale. One can understand this from the fact that the main requirement to 
be satisfied in order the effective theory be approximately SUSY, is that the solution for 
the auxiliary fields of the multiplets being integrated out be always suppressed compared 
with the auxiliary fields of the remaining fields, independently of the dynamics of the 
effective theory. One way of realizing such situation, as found by Brizi et al. [13], is by 
integrating out fields with large SUSY mass compared to the SUSY breaking scale, driven 
this last one by the remaining dynamics. In this case is the large inertia from the SUSY 
preserving sector what shields it from the SUSY breaking effects in the second sector, and 
the auxiliary field VEV's of these multiplets are suppressed by the masses. In our case, is 
the weak coupling among the different sectors what ensures that the solution in one side 
be independent of the other one. Then any back reaction from the SUSY breaking sector is 
negligible in the SUSY preserving sector. As a matter of fact, in case the mixing parameter 
e vanishes the solution for the auxiliary fields exactly vanishes and one can even speak 
about an exact SUSY description of the effective theory, contrary to the case with the 
hierarchy where there are always corrections to the two derivative SUSY description of the 
effective theory [13]. In order to put these arguments in a complete SUSY framework, one 
should understand them from a fully superfield approach, a study left for a forthcoming 
paper [31]. 
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A LVS with matter 



This appendix is devoted to resume the principal features of the LVS models with matter 
and gauge interactions, for details please refer to the original works [21, 29].^^ In the 
following we use the more orthodox form for the scalar potential integrating out the full 
set of auxiliary fields, 

V = e^{GaG^''Ga-3) + ^{Gx?, (A.l) 

with the matrix G"" = G^^, and the rest of the functions as defined in the main text. We 
restrict the study to the simplified version of the models, where the Dilaton and Complex 
structure moduli are regarded as frozen, since our conclusions show that all features we 
need are encoded here. Moreover we neglect the now constant part /Ccs in the Kahler 
potential, eq.(2.15), changing only by a multiplicative factor the whole F-term scalar po- 
tential, and as before we neglect the non-Abelian parts in the L'-term dynamics. 
The Kahler moduli enter in the superpotential only though non-perturbative effects, ex- 
ponentially suppressed by the moduli, then any dependency on the large modulus can be 
neglected. In absence of matter fields the superpotential then looks like 

W = Wo + Ae-''\ (A.2) 

with Wo the 0{1) remnants of a flux induced superpotential once the Dilaton and Complex 
structure are frozen, i.e.. Wo = WsciUo, So), and the amplitude A might also depend on 
these moduli. The standard scenario works without gauge interaction nor matter like 
fields, then the scalar potential, in a large volume expansion and regarding t and all the 
parameters as real, reads [18], 

8aVW-2at \AWo\ , H\Wo\' 

where ^ = ^{2Sr)^^'^. The competing effects of these three terms lead to a stable minimum, 
being crucial that the last term be positive definite, namely > 0, which translates to 
a condition on the kind of CY, in a generic setup, where such minima can exist, i.e., 



^^There is also an interesting realization of LVS in Heterotic compactifications recently developed by L. 
Anguelova and C. Quigley in [51]. However, the phenomenological constrain on the Dilaton VEV, (S) ~ 2 
coming from an good approximate value of the gauge couplings, together with the expression for the string 
coupling gs ~ ^T'^ / S [2], translates to a tough constrain on the VEV of the volume as the theory starts 
now to leave the perturbative regime. 
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as ^ ~ /i^'^ — h^'^ this implies more complex structures than Kahler moduli [23]. The 
minimization of this potential leads to [18] 

with corrections of order 1/at <^ 1. The canonical normalized masses mj- ~ e"^/^ and 
rrit ^ e can be estimated by disregarding the Kahler mixing, i.e., m? ~ {Kii)~^didiV, using 
the notation V^^ ~ e as in the main text. Notice the possibility of realizing low energy 
SUSY breaking, characterized by the gravitino mass = c*^/^ ~ e, with a natural VEV 
for W driven by the exponentially sized volumen. The problem with this setup is the fact 
that the vacuum turns out to be a deep AdS, of ©(e^), as can be seen from the fact that 
the second term in (A. 3) has in the amplitude a factor S> 1 which enhances it compared 
to the other two. 

A.l Charged large modulus 

Once matter fields are considered, with a Kahler potential like the one in eq.(2.18), even if 
these turn out not to appear in the superpotential the F-term scalar potential is affected 
with new terms lead, in a large volume expansion, by 

Vf D e«G««|GQ|2 ~ e^G-^|GQ|2 ~ e'+^\Q\^ . (A.5) 

If such terms are not to sweep off the minimum of the moduli, expected to proceed like 
in the standard scenario, and potentially realize a zero tuned cosmological constant, these 
should scale at most as e^, so we get a constrain on the VEV for the Q's, 

(IQl')<e'-". (A.6) 

A superpotential like eq.(A.2) still is possible with a neutral small modulus, and the ap- 
pearance of the Q fields might be forbidden by gauge invariance, then the generic form for 
the D-term part of the scalar potential is the one give in eq.(3.14). 

Vanishing Q VEV 

The simplest and nicely enough also a realistic scenario developing vanishing cosmological 
constant, is the one where the Q charge has opposite sign of the GS coefficient, so that all 
terms in the L'-term are positive definite, as well the mass term in the F-term potential. 
Then a vanishing VEV solution for Q is possible. 

The uplift, thus, is not possible through (A.5) but from the D-term potential scaling as 



hx \T + TJ hx 

so that Gx ~ and ~ 6'^ e^/'^ /hx- Then in order this to work as an uplift for 

the AdS solution of = —0{c'), and no to spoil the moduli stabilization, we find that 
the GS coefficient should be tuned as 5'^ < e'l^\fhx^ and even in the best case when 
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hx ~ T ~ e~2/'^ we have a strong constrain 6'^ < e^/^. 

Contrary to the standard scenario no massless modes are reahzed, as the imaginary part of 
T is the would-be Goldstone boson eaten up by the now massive vector field of the broken 
U{l)x- For the Q field the mass is ruled by the term (A. 5), so that the canonical mass 
goes like mq ~ e. The gauge breaking scale, encoded in the canonical mass for the vector 

bosons, is given by mx = ■J KrpfX'^ /hx ~ /y/hx ~ e^^^- 



Non- vanishing Q VEV 

When the sign of qq is opposite to the one of the GS coefficient the minimization of the 
potential leads in a good approximation to a cancellation of the D-term potential, driven 
by a non-vanishing Q VEV, so 

im^-^^-^'^'-sr (A.8) 

Let us see the structure for the full e.o.m. for Q, so to have a more precise estimation for 
the VEV of the D-term part of the potential 

OqV « aQ(e«G-^|GQ|2) + -1-GxGqx ~ e^+^Q + -i-e^QGx = 0, (A.9) 

where without loss of generality we regard hx independent of Q, the discarded term being 
of order {Gx)"^ is any way irrelevant. Then, we find that 

Gx ~ hxe^ , r^Gx /hx ~ . (A.IO) 

We have, then, that the L'-term part of the potential scales as Vd ~ hxf^ and even in the 
best case can not work as uplifting for the moduli vacuum which follow again 

like in the pure moduli case. This role is now played by the contribution (A. 5), scaling as 

V.,uft-e^^-\Q\'-e^"5^. (A.ll) 

Again a tuning on the GS coefficient, 5'^ < e^^^, is required though in this case is milder 
and furthermore independent of the scaling of hx- 

The physical mass scalings follow exactly like in the previous case, with an important con- 
tribution to the mass of Q from the D-term dynamics though, and the canonical normalized 
mass for the vector boson mx ~ e^/^/i/^x- 

A. 2 Charged small modulus 

Things are more reach when the small modulus is charged as the the superpotential starts 
necessarily to depend on the Q fields. As is well known depending on the number of families 
and the rank of the gauge group the description of the system changes. 
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Nf < Nc case 

The following model was already proposed by Cremades et al. in [21], however, it was not 
discussed there. Here we discuss the stabilization of the fields in detail and try to argue 
why is not possible to realize Minkowski/dS vacua with this setup. 

Working in the window Nf < Nc for a non-Abelian SU (Nc) sector an ADS nonperturbative 
superpotential [50] is generated, depending on the non-perturbative scale A and the chiral 
fields, Q (Q), transforming in the fundamental (antifundamental) . It is convenient to 
work in D-flat directions characterized by the mesonic fields, then choosing the canonical 
normalized field (f> = \/2QQ [46], and regarding each family as the copy of the others the 
non-perturbative superpotential takes the form 



-Nf 



Wnp = {Nc-Nf)^-^j^j . (A.12) 

The non-perturbative scale is a function of the moduli which as before we take to be the 
small modulus. To make easier the following analysis we take the simplest case, namely 
Nc = 2 and Nf = 1, then the superpotential takes the general form, 

W = Wo- - ^mpcP^ , (A.13) 

where we have included the flux induced superpotential that stabilizes the Dilaton and 
Complex structure moduli, and encoded some Nc and Nf dependencies in the parameters 
A and a that are naturally of order one. A possible mass term is also included, which 
requires a further charged field, p, singlet under the SU{2)}^ 

Without loss of generality we fix the charge for the modulus by (5* = 2/a, so the holomorphic 
Killing vectors are Xx = i(0, 1/a, — <^/2, p). The Kahler potential for the matter fields is 
taken to be if D ~r^\Q\'^^ leading expression for the D-term 

Gx = -^X^G. . i + ^i2p-p - ^4>?j , (A.14) 

where we have also discarded subleading terms in 1/at. 

Despite the mass term for the mesons these cannot be integrated out since the tachionic 
mass from the D-terui dynamics likely dominates. In fact, the minimization of the potential 
leads to a cancellation of the D-term potential via a non-vanishing VEV for cj), analogous 
to the previous case, but now this happens for the opposite sing of the charge as the sign 
of the field dependent FI term changes. Then {(f)'^) w z(^yT^^^'" ~ ^^~^"^^<^o(*f taking 
a priori {(f)) ^ (p) and (poitr) encoding factors scaling like 0(1). In order to stress the fact 
that the dynamics from the mass term are weaker let us discuss first the case m = and 
then see the consequences of turning this term on. Thus, taking the gauge kinetic function 



^^The phase between the terms is irrelevant as can be modified by the value of the coefficient and have 
been fixed only in order to clear up the following formulae meanwhile the coefficients are regarded as 
positive. 
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fx = t the leading D-term potential for p once (p is fixed is 

From here on we consider real solutions for the e.o.m. and forget the subscript r, this is later 
checked numerically for real parameters. Linear terms in p would appear only as factors 
of the mass term with leading expression, coming from the and G^^GtG^ + h.c. 

terms, which once is fixed reads at leading order 

,,p{lE2^^^ + ^aAe-'^'{Ut)-'-atZ\t))) . (A.16) 



X'^Z{t)T3/2 



Therefore, for the case we are currently studying a solution p = is possible. The mass 
terms for p are, at leading order, given by 



2a^e-"*^ {Z'{tf-Z{t)Z"{t)), (A.17) 



again lead by the Gt and terms. With this fixed solutions for p and (j) we get an effective 
potential for the moduli, which at leading order in the volume and 1/at ^ 1 reads, 

_ 8Z(^)2a4A2e-2°* 2yftZ{t)a'^WoAe-''^ 

where the first two terms come from the contribution and the last one is the usual 

a' term from G'^'^\Gt\'^ ■ Interestingly enough the structure of this potential is exactly 
like the standard LVS, eq.(A.3), with a change in the powers for the large modulus. This 
change in the powers induces a different scaling for the exponential since the minimization 
leads to the same scaling of all three terms at the minimum. Indeed we see that at the 
minimum the exponential will in this case scale like e~°* ~ T"~^ ~ g2-2n/3^ while in the 
standard scenario e~"* ~ 2"-3/2 ^ g_ order to have cleaner formulae in the following we 
fix Z{t) = t, then the minimization in T leads to a solution of the form, using at^X, 



^ .27{3 + 2n)Woe''^ I \^ 4(4n - 3)| \ \ \^ ( 27Woe''^ \ 



8(4n-3)a2A 1 \ {3 + 2n)H^/^X ( \A{3 + 2n)Xa^Aty^ , 

(A.19) 

Plugging this solution into the e.o.m. for the small modulus lead by the derivatives of the 
exponential factors one gets the following real solution 

a very similar result to the one of the standard LVS. 

At the true minimun the D-term actually is not zero as the F-term part also contributes 
to the e.o.m for (p. The estimation for its correct VEV is done by refining the e.o.m. for (f> 
including the leading contribution from the F-term potential, analogous to eq.(A.9). This 
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shows that at the vacuum Gx ~ and so it is . 

Like in the standard scenario the vacuum, so far, turns out to be an AdS one, again of 
0{e^)^ since the negative term in (A. 18) leads in a at ^ 1 expansion, and the D-term 
dynamics, although positive definite cannot do the job since scales like 0{e^). 
When a non- vanishing value for the mass parameter is turned on not only a non- vanishing 
VEV for p is generated but also a term of the form, 

^FD^^^^^^mVr (A.21) 

coming from C^IGpp, which being positive definite can potentially uplift the vacuum this 
as far there is tuning in the mass parameter of order m ~ e""^/^. Notice that the term 
(A. 15) might work also since now the VEV of p is non zero. However, the minimization 
of /9, driven by (A. 15) and the linear term leads to that (A. 15), and the 

linear term also, scale like e^^/^ not enough for the uplift. 

Unfortunately by numerical proof we found that the perturbation on the vacuum once the 
mass term is turned on does not allow to uplift and the minimum continue to be an AdS. 
This can be understood from the fact that the uplifting term goes like VupUft ~ 
which cannot compete with the negative term in the moduli potential scaling as \ftc?Z'^e' ^ 
and if one tries to increase the mass parameter in order to compensate this the moduli 
vacuum starts to be perturbed by sending their VEV to larger values so that the net effect 
is lost. We have check numerically, in fact, that the vacuum with a non- vanishing m, but 
still small, ends deeper than in the m = case, then when increasing m only a virtual 
uplift is realized by the shrink of the potential due to a larger T VEV. 
The L'-term dynamics drives the leading contribution to the (\) real part mass, ?7T,^g^^j ~ 
{dfjyGx)"^ ~ <^^/T^'" ~ gi+2n/3^ gQ canonically normalized ~ ^/e^ like the gauge symmetry 
breaking scale. More precisely this real field has a small component from the real part of t, 
and the imaginary counterpart is the would-be Goldstone field eaten up by the now massive 
vector field. The other combination of (j) and t gets a mass like for the small modulus in 
the standard stabilization. For the large modulus we have the same scaling for the mass as 
in the standard scenario. For p we read the mass from the quadratic term eq.(A.17), which 
scales like = g2+2n/3^ g^ normalized rup ~ e for both real and imaginary components. 

Nc <Nf < 3/2iVc case 

An interesting possibility was studied in [29], where standing in the Seiberg duality win- 
dow Nc < Np < 3/2Nc of the non-Abelian sector found vacua with exponentielly sized 
volumina. The dual model is described by a meson like field, (p, and two chiral fields, q 
and p, with a superpotential depending on the non-perturbative scale. A, a possible mass 
term for (j) and a further scale p dictated by the duality relations. As usual the non- 
perturbative scale has an exponential dependency on the small modulus. Then replacing 
the mass parameter by a dynamical field p, we have the following superpotential [29], 

W = Wo + Ae-'''(^^ + p<P^ , (A.22) 

which has an U{1) symmetry characterized be the charges. 
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SO that t gets a non-linear realization identifying (5* = 2/a as the GS coefficient for t. 
With this superpotential and particular modular weights, K D ^|<3p, Krippendorf and 
Quevedo found vacua with exponentially large volumina and vanishing cosmological con- 
stant [29] . In the following we will give a generic study of the scalings around the vacuum. 
To simplify the analysis we neglect the fields q and p, as their VEV turn out to vanish [29] . 
Taking for K D -^fe^lQI^, and fx = t the leading contributions to the D-term potential 

r 

reads ^ 

3 -^+^(2|p|2_|^|2)') , (A.23) 




20^3/2 Ti 



neglecting also terms subleading in 1/at. Like in the previous two cases the minimization 
of the potential would lead to a leading cancellation of the D-term potential due a non- 
vanishing (j) VEV, then \<j)\'^ w ~ e^^'^"'^'^4>o{tr)'^ taking a priori (p) ^ {(j)) and 

'2<iZ {tf )7V 

(j)o{tr) encoding 0{1) factors. Plugging back this solution into the scalar potential one finds 
the following leading potential for regarding real solutions so to forget the r indices, 

(^2^2 -2a i Wo^e-"*, 3CV2 Zitf , ,^ 

« + '^^■^'O + ^ ■ yo«^^^P^. (A^24) 

We have tacitly neglected a term quadratic in p coming from e'^G''^''^|G0|2 D ^n-i^-at ^ 
as can be check at the end, using the resulting scalings, it is indeed irrelevant. Then the 
minimization for p leads to 

^ (2A)2/3 V 2 y Z(t)5/6 T5(3-2n)/12 ' ^^"^^^ 

SO around this solution there is a scalar potential for the moduli from which their stabi- 
lization can be studied, 

^3 ^2g-2at^ 3 9 1/3 (^^2)2/3^2 ^4/3g-4at/3 3Vt/2g 

^2aA2Z(t)2r9/2-2" 2A8/3V2; Z(t)4/3 T^-4n/3 ^ 2X^T^/^ ' ^ ■ ^ 

As was noticed in [29] the structure of the resulting potential is quite similar to the one of 
the standard scenario though the powers of T do not allow for a precise analytical study. 
They find, however, numerically vacua with exponentially large volume which furthermore 
realize vanishing cosmological constant. 

Let us push a bit further the analysis for the stabilization, by writing the potential for the 
moduh as Vmod = (ie-2'^*r2" - ^e-^''*/3j^'^n,/3-i/2 _^ c)/T^/^. By requiring a vanishing 
cosmological constant the e.o.m. for T reads, ^'^ 

drVniod ~ (2nie-2«*r2" - {An/3 - i/2)5e-^'^*/32.4n/3-i/2^ /^ii/2 ^ q ^ (^.27) 



^^The same conclusions are reached using the leading e.o.m. for t, regarding all three factors in the 
potential to be of the same order. 
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so we get that e^"*/'^ ~ j'2n/3+i/2^ telling us that the scaling of the two first terms in the 
moduli potential is the same, more precisely 0{T-^^/^) = 0{e^^/^). But if the cosmological 
constant is to be zero the third term cannot have a larger scaling so one needs to require a 
tuning in the parameters, for example requiring C ~ j"-^/^ e^/^. Such requirement is in 
fact found in [29] by doing an approximate analytic study and looking for conditions for 
the solution to exist. In their approach they consider the possibility of changing a bit the 
powers of T appearing in the potential, to be precise by T^/^, something valid as far the 
volume is not huge. In this regime, then, we can neglect a possible tuning in C or in any 
other parameter appearing in the potential and regard them as 0{1). 
Using the full e.o.m. for (j) we estimate the VEV for the L'-term finding like in the previous 
cases ~ Gx ~ e^- Unfortunately this rough analytic approach is not fine enough so to 
provide a good starting point to look for numerical solutions, then we have been not able 
to check these properties numerically. 

B Numerical tests 

This section shows numerical tests for the decoupling in the models discussed analytically 
in the main text but for the Krippendorf-Quevedo model [29] (see section (A.2)), for which 
unfortunately so far we have not managed to find numerically such vacua. Although the 
form of the chosen Kahler potential and superpotential on each case are string inspired, 
we sacrifice a precise justification of the numerical parameters in favor of avoiding extra 
factors that can make less clear the analysis of the results. 

i?-sector Toy-model 

In order to implement a numerical test we should choose a H sector, which better to 
be as realistic as possible. Our main assumption in the analysis was that the conditions 
dnWsc + WscdnKsc = do not leave flat directions in the H fields. From the way the 
Dilaton enters in the flux induced superpotential of type-IIB superstrings [40], we need 
at least one more dynamical field in order to satisfy such requirement.^^ In the following 
examples we will take the minimal set of one Complex structure field plus the Dilaton 
described by the Kahler potential ICcs = — log(S' + S) — log{U + Lf). The superpotential, 
inspired from explicit calculations on flux compactification [52], is given by 

WcsiU, S) = acsU^ + b,sU + S{fcsU^ + dcsU + ecs) ■ (B.l) 

With numerical parameters that in all following examples we will take as Ocs = 1, bcs = 
dcs = — 3, /cs = 1 and ecs = 5/2, so that lead to the fixed values, solutions to the F- flatness 
conditions. So = Uq = I, and more over Wcs{Uo, So) = Wo = 1. 

The other common feature in the following examples is the matter independent Kahler 
potential, given by the one of type-IIB compactifications with orientifolds and the leading 
a' corrections, eq.(2.15), with the CY volume (2.16). Notice that all these features are 
fairly natural and generic, and no particular intention drives its election. 

Fixing S alone is still possible if one allows a constant part in the superpotential, something that 
obviously we want to avoid. 
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{X) 


A(X) 


px 

^ c 




u 


1 - 4 • 10-20 


-4 • 10-20 


1 ■ 10-41 




s 


1 - 1 • 10-19 


-1 • 10-19 


3 • 10-41 




T 


7.5 • 10^3 


2 • 10-18 


1 • 10-21 


3-10-18 


t 


3.21 


7 • 10-20 


1 • 10-^3 


5-10-18 



Table 1. VEV's of the fields and canonical normalized i^-tcrm, and their relative shifts compared 
to the simplified model. The Q sector is not reported given the triviality of the results. Here and 
in the main text AX = (X — Xsim)/ X. All quantities are in natural units. 

B.l Charged large modulus 

For these models we will take the following particular choice, 

K = K^od + ^jY^^Ynr^\Q\\ W = W,s-AUe-''\ (B.2) 

where we have introduced the mixing between the complex structure moduli and the 
matter fields in the Kahler potential, and a further coupling in the superpotential with 
the small modulus. The Killing vectors for the systems, with ordering {?7, S, T, t, (/>}, are 
{0, 0, z5"^/2, 0, if^gQ}. The gauge kinetic function is chosen iohe fx = T + kS, k a real 
number. 

Vanishing Q VEV 

If the GS coefficient is taken as positive then a negative charge for Q leads to the vanishing 
VEV case presented in appendix A.l and whose decoupling was studied in section 3.2.1. 
We take the following values for the parameters, besides the ones already fixed in the Wcs 
superpotential, 

A = l, A = l, a = 511, i = 2, K = l, 7? = 1, n = \, qq = -1 , 6^ = 1 ■ IQ-" . (B.3) 

The value for the GS coefficient S'^ is chosen such to fine tune the cosmological constant. 
Table 1 shows the solution for the vacuum and the deviation from the simplified model 
where the S and U superfields are fixed to So = Uq = 1, with vanishing spinor and auxiliary 
components. The canonical normalized auxiliary fields VEV, F^ = IGjjF^F^l^^"^ no sum, 
are also shown. The system corresponds to an e parameter, as defined in the main text, 
of roughly O{10~'^^), expecting then corrections to the simplified model of this order. The 
physical mass spectrum in GeV units is 

mt^ = 18.5 - 10-2 , mt^ = 18.4 - 10-2 , mg = 8.8 - 10-2 , mr, = 7.0 - lO'i^ , 
rrig^ = 3.4 - 10-2 , rrig^ = 3.0 - 10-2 , m^^ = 4.5 - 10-^ , m^^ = 8.0 - 10"^ . 

where S and U are defined by a sum of the original 5 and U, but with main component S 
and U respectively, and the subscript r (i) denote the real (imaginary) field components. 
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The imaginary component of T is the would-be Goldstone boson so is exactly massless. 
The shifts compared to the simplified model are 

Amt = 7 ■ 10"^^ , Amq = 5 • 10"^^ , Amr, = 8 • 10"^^ , (B.4) 

defining AX = {X — Xgim)/ X. For the gravitino mass, the D-auxiliary field VEV and the 
cosmological constant we have, all in natural units. 





= 1.3 • 10^^*^ , 


AD"" 


= 5 • 


10-^ 


-18 


"l3/2 


= 7.7- 10^2^ 


Amli^ 


= 3- 


10-^ 


-18 




= -5.5 • 10^^^ , 


AVo 


= 9- 


10-^ 


-18 



Although a quite unrealistic scenario since presents a SUSY breaking scale in the l^~^GeV 
region and a too light moduli, such small e helps to split the scales and clearly shows the 
scaling of the errors in the simplified version. Any way we have check for different values of e 
these scaling. From the results we see that the mistake are few orders of magnitude greater 
than the predicted values, reaching even three orders in worst cases. Still, it is clear that 
all errors are correlated and of the same order. In order to understand such difference with 
the prediction we should remember the numerical factors carefully extracted in section 2.2. 
There we found that the error was enhanced by a factor ~ ^Xt^/'^a^A^ with A a numerical 
factor missed in the analytical solution for e"*. This factor accounts for an increasing of 
two orders in the result. The extra order is easily understood from the fact that the error 
in the small modulus propagates in the e.o.m. for the large one with a factor a. Indeed 
the larger error in the fields VEV corresponds to T, as can be seen in table 1, and all the 
following samples, and then it propagates to the auxiliary fields VEV's, masses and more 
drastically in the cosmological constant, being the first two suppressed by e*^/^ ~ i/7"3/2 
and the last one by ~ l/T^, so that roughly AV ~ —2>V AT. 

Non-vanishing Q VEV 

With the opposite sign for the Q charge the VEV of Q is non- zero and an F-term uplifting 
is possible, section A.l. The parameters taken for this example are 

A = l, A = l, a = 37r, ^ = 2, k=1, r] = I , n=^, qq = 1 , 5^ = ^ ■ 10^^ . (B.6) 

9 

The value for the GS coefficient 5'^ again is chosen such to tune the cosmological constant. 
Table 2 shows the VEV's of the fields and their deviations from the simplified model. The 
system corresponds to an e, as defined in the main text, of roug hly 0(10-^2), expecting 
corrections to the simplified version of this order. The physical mass spectrum in GeV 
units is 

mq^ = 4.2 • 10^ , mt = 3.6 • 10^ , mq^ = 9.8 • 10^ , ttit, = 0.6 

mf^^ = l.l•10^ mJ^^=9.9•10^ m^^ = 1.5 • 10^ m_5^ = 2.6-10^ (B.7) 

again S and U are defined by a sum of the original S and C/, but with main component S 
and U respectively and the subscript r (i) denote the real (imaginary) field components. 
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A / V\ 






u 


1 - 1 • 10-11 


-1-10-11 


1 - 10-24 




s 


1-5 - 10-11 


-5- 10-11 


3 - 10-24 




T 


1.59 • 10^ 


-4-10-10 


4-10-13 


6-10-10 


t 


3.21 


-2- 10-11 


1 - 10-20 


8-10-10 


Q 


1.2 • 10"^ 


6 - 10-11 


1 - 10-20 


8-10-10 



Table 2. VEV's of the fields and canonical normalized _F-term, and their relative shifts compared 
to the simplified model with IS.X = {X — Xsim)/ X. All quantities are in natural units. 

The imaginary component of T is the would-be Goldstone boson so is exactly massless. 
The shifts compared to the simplified model are 

Amt = 5-10-io, Amg^ = 4 - 10-10 , Amg^ = 6 - lO-io , Amr, = 6 - lO-io . (B.8) 

For the gravitino mass, the Z)-auxiliary field VEV and the cosmological constant, in natm'al 
units, we have 





= 2.1 - 10-26 , 


AD^ = 1 


10-^ 




"i3/2 


= 2.5 - 10-13 , 


Am2/2 = 6 


10-10, 


(B.9) 


^0 


= -4.3 - 10-39 , 


Ayo = 1 


10-^ 





This, now more realistic, scenario, with a SUSY breaking scale of 0{lQ^)GeV presents 
exactly the same features in the errors as the previous case. This is expected as the 
largest corrections come from the same source independent of the matter fields and gauge 
dynamics. Thus, the explanation for the factors showing the errors larger than expected 
ones follows verbatim here. 

For the record we have done numerically the very same model changing the gauge kinetic 
function to / = t + kS. This in order to check the claim that is the scaling of the VEV 
for auxiliary field what matters. Otherwise, as explained in footnote 10, with f ^ t 
the corrections would rather scale like 0{e^/'^). We find in fact numerically the very same 
results as the ones just presented, being the only drastic change the VEV of Gx- 

B.2 Charged small modulus 

The system is the one described in detail in appendix (A. 2), so besides checking explicitly 
in numbers the decoupling we take the chance to see all features found analytically there. 
For the matter fields the Kahler potential is taken to be 

KQ = jr{\p? + Ur\4>?) . (B.IO) 

The superpotential is taken as 

W = Wcs{U,S)- AS^-mU^ p4>'' , (B.ll) 
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1 - 4 • 10-5 


-4- 10-5 
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12 
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11247 


-2- 10-^ 


7 


• 10' 


-7 


3 • 10-^ 
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2.36 


-2- 10-5 


4 


10- 


11 


4 • 10-4 




4.6 • 10-2 


7 • 10-5 


9 


10- 


11 


4 • 10-4 


P 


7.1 • 10-3 


8 • 10-5 


2 


10- 


10 


3 • 10-4 



Table 3. VEV's of the fields and canonical normalized F-terms, and their relative shifts, derived 
by a numerical analysis with AX = {X — Xsim)/X. All quantities are in natural units. 

and the gauge kinetic function fx = t + i^S, which is different from the one taken in (A. 2) 
but its t dependency and scaling in e are the same though. The parameters taken are 

6 2 1 

A = l, A = l,a = 37r,^=-,K = 1, m = 268 • 10-5 , 5' = — = -- ,qp = I . 

(B.12) 

We take a value for m that in principle would work to uplift the vacuum, being 0(ei/2). 
However the relative uplift to the zero m vacuum turn to be only a virtual one being only 
due to the shift caused on the VEV of T which moves to larger values. In fact, as advertised 
before, for m slightly smaller than e^^"^, such that its perturbation on T is less relevant the 
vacuum goes slightly deeper than in the zero m case, the same happening with positive or 
negative values for m. As defined in appendix (A. 2) the GS coefficient is related to the 
parameter a by 5* = 2/a, and the (p and p charges are fixed as well, but we report them 
again for completeness. 

Table 3 shows VEV's of the fields and the deviation from the simplified model. The 
canonical normalized auxiliary fields VEV are also shown. The system corresponds to 
an e, as defined in the main text, of roughly 0(10-^), to be precise 8 • lO-'', expecting 
corrections to the simplified model of order ©(e^/^) 10-5. We check then that indeed the 
field (j) is stabilized mainly by the nearly vanishing D-term condition, (j) ~ (a^t^T,.)-^/^^ pgr 
the rest of the quantities the scalings are also in agreement with the one found in appendix 
(A. 2). The physical mass spectrum in GeV is 

m^^ = 6.3 • 10^4 ^ ^, ^ 4 5 . ^qIS ^ ^ 2.2 • 10^^ ^ ^ 2.O • 10^^ ^ 

rufj^ = 1.9- 10^^ , m^^ = 1.7 -10^3^ m^, = 2.5 • 10^^^ m^^ = 4.4 • 10^\ (B.13) 

mr, = 1.2 • lO'^ . 

where (p (t) is a linear combination of (p and t with (j) (t) as main component, and S and U are 
defined by a sum of the original S and U, but with main component S and U respectively. 
The subscript r (i) denote real (imaginary) field components. The imaginary component 
of 4> is the would-be Goldstone boson so is exactly massless, as well the imaginary part of 
T. 

The shifts compared to the simplified model are 

Am, =2-10-4, Amj = 3-10-4, Am^ = 2 • 10-4 , Amr, = 3 • 10-4 . (B.14) 
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For the gravitino mass, the D-auxihary field VEV and the cosmological constant we have, 
in natural units, 

= 4.5 • 10"^^ , AD^ = 5 • 10"^ , 

my2 = 4.2 • 10"^ , ^^1/2 = 3 • lO""^ , (B.15) 

Fo = -1.7- 10-20, AVo = 7-W-'^. 

This scenario with high SUSY breaking scale, 0(10^^)06^, shows an interesting situation 
that might happen with the corrections. Although the corrections fall in predicted order, 
i.e., ©(e^/^), the main corrections come from the 0{e) contribution coming again mainly 
from the G(f|F*p terms. This is check by performing the same numerical search but 
setting m = 0. Analytically it can be understood by looking to enhancement factors 
as we did before to explain the numerical results for the previous two cases. For the 
corrections coming from the mass parameter we have A ~ " VjAm ^2/3^ with A = e~"'^ /e^^^, 
m = m/e^/2, p = p/e^^'^, and for the present example the correction factor is ~ 1/2. On 
the other hand, the ones from Gu ~ ae~"'^/(j)'^ + e^^^4>'^ leads to corrections going like 
A ~ a^t^/"^ A^e, which enhance them by almost two orders. The result is that both errors 
enter in the same order though the ones scaling as 0{e) turn out to be slightly bigger. 
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